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Предложен метод параллельной генерации таких комбинаторных объектов, как перестановки и сочетания. В основе метода ле- 
жит возможность нумерации перечисляемых объектов так, что по номеру можно легко восстановить сопоставленный ему объект. 
Приведены формулы вычисления номера заданного объекта (перестановки, сочетания) и алгоритмы построения объекта по его 
номеру. Представлены результаты экспериментов, подтверждающие эффективность метода. 


В проектировании управляющих систем возника- 
ет необходимость в генерации таких комбинаторных 
объектов, как перестановки и сочетания. Например, 
в задачах проектирования дискретных схем на базе 
программируемых больших интегральных схем с мат- 
ричной структурой возникает задача генерирования 
всех допустимых в некотором смысле подмножеств 
заданного множества. Их генерацию можно вести па- 
раллельно на вычислительных узлах многопроцес- 
сорной системы. Для этого требуется задать нумера- 
цию перечисляемых объектов так, что по номеру 
можно легко восстановить сопоставленный ему объ- 
ект. Построенная таким образом нумерация позволя- 
ет разбить все множество объектов на классы, мощ- 
ности которых отвечают производительности про- 
цессоров в системе, и генерировать объекты различ- 
ных классов параллельно на соответствующих про- 
цессорах. Предлагаются формулы вычисления номе- 
ра заданного объекта (перестановки, сочетания) и ал- 
горитмы построения объекта по его номеру. 

Генерация перестановок 

Рассматриваются перестановки (а„...,а„) из эле- 
ментов множества ІѴ={1,2 т.е. перестановки 
степени я. Известен алгоритм генерации (перечис- 
ления) всех перестановок в лексикографическом 
порядке [1], т.е. так, что перестановка (а ь ...,а„) 
предшествует перестановке (Ь и Ь ъ ...,Ь п ), если для не- 
которого і, 1</<я, имеет место а х =Ъ { , а 1 =Ь 1 ,...,а=Ь. и 
а іА <Ь іА . Генерация выполняется путем преобразова- 
ния текущей перестановки к следующей. Напри- 
мер, перестановки степени 3 перечисляются в сле- 
дующем порядке: 123, 132,213,231,312, 321. 

Назовем индексом перестановки ее порядковый 
номер / в последовательности всех перестановок 
степени л, расположенных в лексикографическом 
порядке, 1</<л!. Рассматриваются задача вычисле- 
ния индекса заданной перестановки и обратная ей - 
задача построения перестановки по ее индексу. Пер- 
вую можно решить путем перечисления и последо- 
вательной нумерации перестановок до тех пор, пока 
не будет получена заданная перестановка, что край- 
не неэффективно. Предлагается вычислять индекс / 
перестановки (а ь . по следующей формуле 

1 = (X («' - ОО - О 0 + (а'п -і - (и - 2)), 

1=1 


где а'=а,+& и к - есть количество элементов а і та- 
ких, что о ( >о, и /</. 

Пример. Пусть я=4, (а 1 ,...,д 4 )=(2,3,1,4). Тогда: 
(й[,й 2 ',й 3 ',й 4 ')=(2, 3,3,4), /=(2-1)3!+(3-2)2!+(3-2)=9. 

Покажем каким образом перестановку можно 
построить по ее индексу. Пусть 0<«п и А=(а ь ..., а,) 
есть префикс перестановки с индексом 1. Обозна- 
чим число перестановок, предшествующих всем 
перестановкам с префиксом А, и положим /V/, =/-/.. 
Назовем М, кратностью префикса А. По определе- 
нию М,~І и /V/, есть число перестановок с префик- 
сом А, индексы которых не превосходят /. Пусть 
У={Уь-,Уп-м}=М-{аі,-,а,-і} и у, Тогда по 

свойству перестановки а,еУ. Допустим, что д =у ь 
\<к<п-1+\. Тогда А=(а ь ...,а 1А ,у к ) и /,=/ м +(/с-])(я-1)!, 
поэтому М=І-ф=І-У, -(к-\)(п-\у~М, -(к-\)(п-()\. 
Таким образом, зная индекс перестановки и любой 
ее префикс, можно вычислить кратность последнего. 

Теорема 1. Пусть 1<Кл-1, (д ь . ..,а, А ) - префикс 
перестановки с индексом I и М іА - его кратность. 
Пусть также {у 1 ,...,у п _„ і }=П-{а 1 ,...,а,-і} и у,<...<у лЧ+1 . 
Тогда /-Й элемент а, в перестановке с индексом / 
есть у к для 


к = 


М,_, + (п — 1)\— 1 
(п-І)І 


Доказательство. Пусть С. для \<к<п—1+\ есть 
множество всех перестановок с префиксом 
А=(а 1 ,...,а 1А ,у к ). Число перестановок в нем |С*|=(л-/)!. 
Ввиду у к <у к+] все перестановки в С к предшествуют 
всем перестановкам в С м . Множества С Ь С 2 , ...,С„_ (+1 
образуют разбиение множества всех перестановок с 
префиксом (а ь ...,а, А ), и перестановка с индексом / 
входит в одно из них. Пусть / обозначает число всех 
перестановок, предшествующих перестановкам с 
этим префиксом. Тогда /Ц ,=/-/ и индексы всех пе- 
рестановок в С к суть /+(&-1)(я-?)!+1, Ф+(к- 1)х 
х(л-?)!+2,...,/+А:(я-0!+1. При (А:-1)(я-?)!+1< 
<М іА <к(п-і)\ перестановка с индексом /входит в С к 
и а,=у к . Прибавив ко всем частям последнего нера- 
венства величину (л— /)!— 1, получим (к-1)(п-1)\+ 
+{п-і)\<М ІЛ +{п-і)\-\<к(п-і)\+(п-і)\-\, или к{п-і)\< 
<М 1А +(п-1)1-К(к-1)(п-1)\, откуда 


к ^ М, А + (п-і)\-\ 
(»- 0 ! 


< Ач-1. 
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Естественные науки 


Следовательно, к 


М і] + (п - 1 ) !- 1 

( п-г)\ 


Теорема 1 доказана. 

При известном префиксе (а ь ...,а„ А ) перестанов- 
ки с индексом / ее последний член находится как 
единственный элемент множества /Ѵ-{а,,...,а„_]}. 

Пример. Пусть п= 4. Построим перестановку 
(а ь а ъ а ь а 4 ) по индексу / =9. 

1=1: М 0 =9, У={ 1, 2, 3, 4}, к= 2; следовательно, а,=2. 
1=2: М,=М 0 -(Л-1)(л-1)!=3; У={ 1, 3, 4}, к= 2; следо- 
вательно, а,=3. 

1=3: М г =М 1 — (А:— 1)(л— 2)!=1; Т={1, 4}, &=1; следова- 
тельно, д 3 =1. 

1=4: Т={4}; следовательно, а 4 =4. 

В результате получаем перестановку (2, 3, 1,4). 

Данный метод вычисления перестановки 
(а ь ...,а„) по ее индексу / можно выразить следую- 
щим алгоритмом. 


Алгоритм 1 (вычисление перестановки (а ь ... 
по индексу Г) 

1. 1=1;(у 1 у 1 ...у„)=(12...п). 


2 . 


Если 1=«, то п. 5; иначе 


/ + (77 - Г) !- 1 

(п-і)\ 


и 


А,) 


3. а,=Уь ІУѵУп-і)=(УѵУк-іУт-Уп^іУ, 

4. І=І-(к-\)(п-і)\\ і=1 + 1 перейти к п. 2. 

5. а=уі. 

Реализация алгоритма 1 предполагает исполь- 
зование "больших" чисел (целых чисел, двоичное 
представление которых требует большего числа 
разрядов, чем в стандартных типах) для задания 
индекса. Вместе с тем, опираясь на него, можно 
предложить следующий алгоритм параллельной ге- 
нерации перестановок. 

Пусть требуется построить все перестановки сте- 
пени п на однородной многопроцессорной вычисли- 
тельной системе, состоящей из р одинаковых про- 
цессоров с общей или распределенной памятью. 
Максимальную загруженность процессоров можно 
обеспечить, разделив последовательность всех пе- 
рестановок на блоки по $=[п\/р\ элементов в каждом. 
Блок с номером т начинается с перестановки, соот- 
ветствующей индексу 1=${т-\)+\ , и генерация пе- 
рестановок в нем возлагается на /и-й процессор. Ес- 
ли число всех перестановок не кратно р, то процес- 
сор, генерирующий перестановки последнего блока, 
строит также (и! той р) оставшихся перестановок, 
число которых меньше р. При выполнении этого ал- 
горитма отсутствует взаимодействие между процес- 
сорами, и все процессоры одинаково загружены, что 
позволяет добиться ускорения (величина, показыва- 
ющая, во сколько раз параллельный алгоритм рабо- 
тает быстрее последовательного) близкого к р. 

Описанный параллельный алгоритм обладает 
высокой эффективностью благодаря пропорцио- 
нальному распределению вычислений и отсут- 
ствию обменов данными между процессорами. 


Предложенный алгоритм был реализован с по- 
мощью технологии параллельного программирова- 
ния МРІ (Ме$8а§е Ра$ып§ ІШегІасе), предназначен- 
ной для организации параллельных вычислений в 
системах с распределенной памятью, и исследован 
в эксперименте на однородной многопроцессор- 
ной системе кластерного типа. На рис. 1 представ- 
лены результаты вычислительного эксперимента, 
из которых следует, что средняя эффективность ис- 
пользования системы алгоритмом равна 0,94. 



Рис. 1. Ускорение для алгоритма параллельной генерации 
перестановок степени п, достигаемое при использо- 
вании 2, 3, ..., 1 6 процессоров 

Генерация сочетаний 

Здесь под сочетанием понимается всякое упоря- 
доченное по возрастанию к- элементное подмноже- 
ство «-элементного множества N={1, 2, ..., п}, т.е. 
любой такой вектор (Ь ь Ъ ъ ..., Ъ к ) длины к<п с ком- 
понентами в И, называемый также сочетанием из п 
по к, в котором Ъ,<Ь І+ 1 для 1= 1, ..., к- 1. В [2] описан 
алгоритм генерации сочетаний, в котором сочета- 
ния порождаются в лексикографическом порядке, 
т.е. так, что сочетание (Ь ь Ъ ъ ..., Ъ к ) предшествует со- 
четанию (а,, а 2 , ..., а к ), если для некоторого 1, 1 <і<к, 
имеет место Ь { =а ъ Ь 2 =а ъ ..., Ь =а, и Ь м <а м . Напри- 
мер, все сочетания из 5 по 3 перечисляются в следу- 
ющем порядке: (1, 2, 3), (1, 2, 4), (1, 2, 5), (1, 3, 4), 
(1, 3, 5), (1, 4, 5), (2, 3, 4), (2, 3, 5), (2, 4, 5), (3, 4, 5). 

Индексом сочетания назовем порядковый номер 
/этого сочетания в последовательности всех соче- 
таний, расположенных в лексикографическом по- 
рядке, 1 </<С*. Индекс /сочетания ( Ь ь Ъ ъ ..., Ь к ) из п 
по к может быть вычислен по следующей формуле: 

/ = ! + Е X Ся г де Ъ» = 0 и сѵ = I . 

і=1 і=Ьі -\+ 1 

Пример. Вычислим индекс сочетания (2, 4, 5) из 
5 по 3. 

1 = 1+ X С 4 2 + ХС‘+ X С 0 ° = 1+ 6+ 2+ 0= 9. 

у'=0+1 у=2+1 ./'=4+1 

Рассмотрим теперь, как сочетание (Ьі,Ь 2 ,...,Ь к ) 
строится по его номеру /. Пусть \<1<к и 
0 =Ь 0 <Ь 1 <Ь 1 <...<Ь,_ 1 . Для любого / в {Ъ, А +\,Ъ, А +2,..., 
п-к+Г}, через обозначим число всех сочетаний с 
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префиксами (Ъ ь Ъ ъ ...,Ъ іА ,г), где ге{Ь )А +1,Ь ІА +2,..^'-Ц. 
В каждом из этих сочетаний недостающие к-1 эле- 
ментов берутся из п-г элементов гМ, поэтому 



Г = Ь ,_! + 1 *=1 


ло всех сочетаний с префиксами (Ъ ь Ъ ъ ...,Ъ ІА ,а) для 
і <(- 1 и а,<й„ то есть всех сочетаний, предшествую- 
щих сочетаниям с префиксом (Ь ь Ь 2 ,...,Ъ 1А ). 

Теорема 2. Пусть 1 <1<к и {Ь 2 ,Ь 2 ,...,Ь, ,) есть пре- 
фикс сочетания с индексом /; тогда 1-й элемент Ъ, 
этого сочетания равен наименьшему из таких у, что 
Ь іА <}<п-кМ и если такие у существуют, и 

Ь,=п-к+і в противном случае. 

Доказательство. Для любого у, где Ъ 1А <]<п-к+1, 
введем І=к+а]. По определению 5 1 и 5" есть число 
всех сочетаний, предшествующих сочетаниям с 
префиксом (Ъ ь Ъ ъ ...,Ь, А У), т.е. /) является индексом 
последнего сочетания с префиксом (Ъ ь Ъ ъ ...,Ъ, А ^- 1), 
а І т - индексом последнего сочетания с префик- 
сом (Ь ь Ъ ъ ...,Ь, А У). Рассмотрим возможные случаи, 
имея в виду, что Ь ІА <Ъ,<п-к+і. 

1. Не существует такого у, что Ъ 1А <]<п-к+1. Тогда 
Ъ=п-к+і. 

2. Для всех у, где Ь іА <]<п-к+1, имеет место I >І ]А . В 
этом случае Ъ,>] для всякого у в границах 
Ъ, А <]<п-к+1 и потому Ь,=п-к+і. 

3. Существуют у с ограничениями Ъ, А <і<п-к+і и 
1<1 }А . Для любого из них /)</. Возьмем наименьшее 
/ из таких у. Тогда І< I и Ьу/. Следовательно, Ъ=]’ . 

Теорема 2 доказана. 

Пример. Построим сочетание из 5 по 3 с индек- 
сом 7=9. 

(=1:х'=0; У'=1, 5' +І =4=6, 5'+5' +1 =6<7; у=2, 4=9, 

У+4+і = 9>7; следовательно, 6,=2. 

/=2: У=4=6; у-3, 4 +і = 4 = 2, ■?' +4 + ,=8</; у=4=и-А:+/; 
следовательно, 6 2 =4. 

/=3: не существует у для 4=Ь, А <]<п-к+1=5, поэтому 
Ь 2 =п-к+і= 5. 

В результате получаем сочетание (2, 4, 5). 
Данный метод вычисления сочетания из п по к 
по его индексу /можно выразить следующим алго- 
ритмом. 

Алгоритм 2 

1. Ь 0 = 0, 5=0, 1=1. 

2. Если 1<к, то п. 3, иначе п. 8. 

3- І=Ь,- 1+1. 

4. Если у <п-к+1 и .ѵ+С„4'< /, то п. 5, иначе п. 6. 

5. 5=5+С*-';у=у+1; п. 4. 

6 . Ь,= % 

7. ?=?+ 1 ; п. 2. 

8. (4, ■■■, Ъ к ) - искомое сочетание. 
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Параллельная генерация сочетаний осущес- 
твляется так же, как и перестановок, а именно: все 
множество сочетаний делится на блоки, мощности 
которых отвечают производительности процессо- 
ров в системе; в каждом блоке по алгоритму 2 с по- 
мощью индекса вычисляется первое сочетание, на- 
чиная с которого один из процессоров системы ге- 
нерирует сочетания всего блока. Как и в случае пе- 
рестановок, благодаря пропорциональному расп- 
ределению вычислений и отсутствию обменов дан- 
ными между процессорами, параллельный алго- 
ритм обладает высокой эффективностью. 

Результаты эксперимента представлены на рис. 2, 
где изображен график усредненного значения ускоре- 
ния, полученного в примерах со следующими парамет- 
рами: и=30, ^=15; я=40, &=10; я=40, к=\5; л=40, &=30; 
и=50, к= 10; л=50, ^=40; й=100, к= 6; я=100, к=94. Сред- 
няя эффективность использования системы (отноше- 
ние ускорения к числу процессов) равна 0,87. 



Рис. 2. Среднее ускорение для алгоритма параллельной ге- 
нерации сочетаний, достигаемое при использовании 
2, 3, 1 6 процессоров 

Отметим, что в случае генерации сочетаний эф- 
фективность распараллеливания меньше, чем в 
случае генерации перестановок. Причина этого в 
том, что время, требуемое на построение блока со- 
четаний, зависит от того, какие именно сочетания 
в нем генерируются. В ходе эксперимента также 
было установлено, что процессоры, генерирующие 
сочетания первых блоков, заканчивают работу 
немного быстрее, чем процессоры, генерирующие 
сочетания последних блоков. Таким образом, при 
генерации сочетаний с большими индексами тре- 
буется, как правило, больше операций, и поэтому 
блоки одинакового размера могут строиться за раз- 
личное время, однако для практических приложе- 
ний эта разница не существенна. 

Полученные в эксперименте значения коэффи- 
циентов эффективности 0,94 и 0,87, соответственно 
для алгоритмов генерации перестановок и сочетаний, 
подтверждают хорошую масштабируемость прог- 
рамм, то есть сохранение тенденции уменьшения 
времени работы при увеличении числа процессоров. 


2. Липский В. Комбинаторика для программистов. — М.: Мир, 
1988.-200 с. 
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